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We study the Dalitz plot of the decay D + — > K~n + n + with a sample of 15090 events from 
Fermilab experiment E791. Modeling the decay amplitude as the coherent sum of known Kn 
resonances and a uniform nonresonant term, we do not obtain an acceptable fit. If we allow the 
mass and width of the A"q(1430) to float, we obtain values consistent with those from PDG but 
the x 2 P er degree of freedom of the fit is still unsatisfactory. A good fit is found when we allow 
for the presence of an additional scalar resonance, with mass 797 ± 19 ± 43 MeV/c 2 and width 
410 ± 43 ± 87 MeV/c 2 . The mass and width of the A5(1430) become 1459 ± 7 ± 5 MeV/c 2 and 
175 ± 12 ± 12 MeV/c 2 , respectively. Our results provide new information on the scalar sector in 
hadron spectroscopy. 

PACS numbers: 13.25.Ft 14.40.Ev 



In this paper we present a Dalitz plot analysis of the 
Cabibbo-favored decay D + — > A'~7r + 7r + using data from 
Fermilab experiment E791. Previous analyses of this de- 
cay jj| |^] modeled the amplitude as the coherent sum 
of known Kir resonances and a uniform non-resonant 
(NR) term. They observed that the NR term is strongly 



dominant, unlike other D decays, and that the sum of 
the decay fractions substantially exceeds unity, indicat- 
ing large interference. Moreover, the fits did not describe 
the Dalitz plot distributions well. In our analysis, we ob- 
tain similar results but with higher statistics. We study 
variations in the underlying model, including changes in 
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form factors, timing of resonance parameters, and the 
addition of known and new resonance structures. In par- 
ticular, we investigate the scalar sector for which there 
has been much uncertainty for many years. 

This study is based on the Fermilab E791 sample of 
2 x 10 10 events produced from interactions of a 500 GeV/c 
tt~ beam with five thin target foils (one platinum, four 
diamond). Descriptions of the detector, data set, recon- 
struction, and vertex resolutions can be found in Ref. @. 
A clean sample of K~tt + tt + decays (charge-conjugate 
modes are implicit throughout this paper) was selected 
by requiring that the 3-prong decay (secondary) vertex 
be well-separated from the production (primary) vertex 
and located outside any solid material. The sum of the 
momentum vectors of the three tracks from the secondary 
vertex was required to point to the primary vertex, and 
each of the three tracks was required to pass closer to 
the secondary vertex than to the primary. We restricted 
the "p\ and xp ranges of the D + candidates to ensure an 
accurate model of our experiment in the Monte Carlo 
(MC) simulation. Finally, we required that the odd- 
charge track (track with charge opposite that of the 
candidate) from the secondary vertex be consistent with 
kaon identification in the Cerenkov counters Q. 

We fit the K~tt + h + invariant mass distribution shown 
in Fig. 0(a) by the sum of D + signal and background 
terms. The signal was represented by the sum of two 
Gaussians, with parameters determined by the fit. We 
used MC simulations and data to determine both the 
shape and the size of charm backgrounds. The signifi- 
cant sources are reflections from Df — > K~K + tt + (via 
K*K + and (j>TT + ), in which one kaon is misidentified as a 
pion. Other sources of charm background are either neg- 
ligible or broadly distributed and thus safely included 
when we estimate combinatorial background, which was 
represented by an exponential function. The number of 
D + candidates obtained from the fit is 16190 ± 139. 

For the Dalitz plot analysis, we selected candidates in 
the Ktttt mass range 1.85-1.89 GeV/c 2 (crosshatched re- 
gion m Fig. 0(a)). This results in 15090 events, with 
about 6% due to background. Fig. 0(b) shows the corre- 
sponding Dalitz plot, m 2 2 vs. m 2 3 , in which the kaon can- 
didate is labeled particle 1, and the plot is symmetrized 
with respect to the two pions (particles 2, 3). 

To study the resonant structure in Fig. 0(b), an un- 
binned maximum likelihood fit is used. The likelihood 
C is computed as C = ILvents 2i=i n BiV Bi + n s Vs 
where and Vs are the normalized probability den- 
sity functions (PDF's) for background and signal, re- 
spectively, and riBi and ns are their fractional contri- 
butions. Each background PDF is written as V Bi = 
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FIG. 1: (a) The K~tt + it + invariant mass spectrum. The 
shaded area corresponds to the background level. The 
crosshatched region is the sample used for the Dalitz plot 
analysis, (b) The D + — > K~tt + tt + Dalitz plot, symmetrized 
for the indistinguishable pions. 



to events above the signal peak in the Ktttt mass range 
1.92-1.96 GeV/c 2 . The shapes of the D+ -> K*K+ and 
Df — > <j)TT + backgrounds are from MC simulations. 

The signal PDF is V s = ^g(M)e(m 2 12l m 2 13 )\A\ 2 , 
where N$ is the normalization, g(M) describes the sig- 



nal shape in the Ktttt mass spectrum, and e(m 
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the acceptance across the Dalitz plot, including smearing. 
The signal amplitude A is a coherent sum of a uniform 
NR amplitude and resonant Kit amplitudes, 
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where each term is Bose symmetrized for the pions: 
An = -4„[(12)3] + .An [(13)2]. The coefficients a n are 
magnitudes and the S n are relative phases. 

Our first fit, referred to as Model A, includes only 
well-established resonances and fixes their masses and 
widths to PDG || values. This approach has been used 
in previous Dalitz-plot analyses (e.g., Refs. g, [|). The 
NR amplitude Ao is represented by a constant; i.e., it 
has no magnitude or phase variation across the Dalitz 
plot. Each resonant amplitude A n {n > 0) is writ- 
ten as A n = BW n F ( D T) F n J) M [ n ] . The BW n fac- 
tor is the relativistic Breit-Wigner propagator, BW n — 
j^-b i (M)J r B i (ml 2 ,ml 3 ), where N Bi is the normaliza- {m 2 — m 2 — im„r„(m)} _1 , where m is the invariant 

mass of the Ktt pair forming a resonance (either m\ 2 
or TO13), m n is the resonance mass, and T n (m) is the 



tion, bi(M) is the distribution in the Ktttt mass spec- 
trum, and J^Bi is the shape in the Dalitz plot. The shape 
of the combinatorial background is obtained from a fit 



mass-dependent width. The factors and F r \ J> are 
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Blatt-Weisskopf penetration factors |7j] , which depend on 
the spin J and the radii of the relevant mesons. In Model 
A, the radii are fixed as ro = 5 GcV -1 for the D meson 
and rji — 1.5 GeV -1 for all Kn resonances 0. No form 

factors F are used for scalar resonances. The term M ( n J) 
accounts for the decay angular distribution. Ref. Q gives 
detailed expressions for all these functions; note that we 
use the opposite sign for the BW n term, for easier com- 
parison of our results to those of Ref. Q . 

For Model A, we fix the NR parameters to be ao = 1 
and Sq = 0, and include all well-established Kir reso- 
nances; the only free parameters of the fit are the magni- 
tudes a n and phases S n of the resonances. The so-called 
decay fraction for each mode is obtained by integrat- 
ing its intensity (squared amplitude) over the Dalitz plot 
and dividing by the integrated intensity with all modes 
present. The fit results are listed in Table [| We observe 
contributions from the same channels reported previously 
[0, |^]; i.e., a high NR decay fraction (over 90%), followed 
by K*(U30)n+, K*(892)tt+, and if*(1680)vr+. We also 
measure a small but statistically significant contribution 
from K$ (1430)7r + . No other resonances considered are 
found to contribute. The sum of the decay fractions is 
~ 140%, indicating a high level of interference. 

To assess the quality of the fit, we developed a fast-MC 
algorithm which produces binned Dalitz plot densities 
according to signal and background PDF's, including de- 
tector efficiency and resolution. A x 2 is calculated from 
the difference between the binned Dalitz-plot-density dis- 
tribution for data and that for fast-MC events generated 
using the parameters obtained from the fit of Model A. 
The x 2 summed over all bins is 167 for 63 degrees of 
freedom (v). The largest contributions to this x 2 come 
from bins at low Kit mass. In Fig. ||(a) we show the 
mass-squared projections; the top (bottom) plot shows 
the lower (higher) mass combination. The points rep- 
resent data and the solid line represents fast-MC simu- 
lation of Model A. The main discrepancies occur below 
0.6 (GeV/c 2 ) 2 and around 2.5 (GeV/c 2 ) 2 . These discrep- 
ancies, and the large value of x 2 l v i motivated us to study 
alternative ways to model the decay amplitude. 

For our second fit, Model B, we allow the mass and 
width of the scalar Aq(1430) resonance to float. In addi- 
tion, we include form factors to account for the finite size 
of the decaying mesons in this scalar transition || [To| . 
The amplitude is written as Fn BW n , in which the 
form factors are Gaussian: = exp(-p* 2 /(2fc 2 )). 

The factor p* is the momentum of the decay products, 
ko = \/6/r, and r is the decaying meson radius. These 
radii (ro and introduced above) become additional 
free parameters in the fit. The results of this fit are listed 
in the middle column of Table |. The decay fractions ob- 
tained are very similar to those found for Model A, but 
the x 2 l v is improved, dropping from 167/63 to 126/63. 
The mass and width of the Kq (1430) obtained by the fit 



TABLE I: Results of the Dalitz plot fits. Models A and B 
are without k; Model C is with k. For each mode the first 
row lists the decay fraction in percent, the second row lists 
the magnitude of the amplitude (a n ), and the third row lists 
the relative phase (S n ). The first error listed is statistical, and 
the second error (when listed) is systematic. 
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are 1416 ±27 MeV/c 2 and 250±21 MeV/c 2 respectively, 
which are close to the PDG values of 1412 ± 6 MeV/c 2 
and 294 ± 23 MeV/c 2 ||. The meson radii obtained are 
r D = 0.8 ± 1.0 GeV" 1 and r R = 1.8 ± 3.4 GeV" 1 . 

Since Model B still does not give a satisfactory fit, we 
allow for an additional scalar amplitude (Model C). For 
this extra amplitude, we use Gaussian form factors sim- 
ilar to those used for the A* (1430) [JllJ . The decay frac- 
tions and relative phases obtained by the fit are listed in 
the right-most column of Table |. In the table we denote 
the additional scalar resonance as "k". In fact, discus- 
sions of the existence of such a resonance are found in 
the literature jl2], [L3|. The fit results are very different 
from those obtained for Models A and B, the NR decay 
fraction drops from 90% to (13 ± 6)%; the ktt + channel 
is dominant with a decay fraction of (48 ± 12); and the 
sum of all fractions is ~ 90%, with smaller interference 
effects. The x 2 l v decreases to 46/63, substantially lower 
than those for Models A and B. For the A* (1430) res- 
onance, we measure m^*(i4 3 o) = 1459 ± 7 ± 12 MeV/c 2 
and r K - (1430) = 175 ± 12 ± 12 MeV/c 2 . These val- 
ues are significantly higher and narrower, respectively, 
than those given by the PDG Q which are taken from 
LASS @. See also @. The mass and width of the 
additional resonance (k) are 797 ± 19 ± 43 MeV/c 2 and 
410 ± 43 ± 87 MeV/c 2 , respectively. The meson radii 
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FIG. 2: m 2 (Kir)i ow and m 2 {K-k)^^ projections for data 
(error bars) and fast MC (solid line): Models A (a) and C 
(b). 



obtained in Model C are ro = 5.0 ± 0.5 GeV -1 and 
Tr = 1.6 ± 1.3 GcV -1 . The Kit mass-squared projec- 
tions are shown in Fig. ||(b). 

To better understand our results for Model C, we per- 
form the following test. For Models B and C we use the 
fast-MC to generate an ensemble of 1000 "experiments," 
with each experiment having a sample size Poisson- 
distributed around our observed sample size. For each 
experiment we calculate Awb,c = — 2(ln£s — hi£c), 
where Lb and Cc are the likelihood functions evaluated 
with parameters from Models B and C, respectively. For 
the ensemble generated without ktt + , (Awb,c) = —123; 
i.e., Model B has greater likelihood. For the ensemble 
generated with ktt + , (Awb.c) = 143; i.e., Model C has 
greater likelihood. In both cases the rms of the distri- 
butions is about 23. For the data, Awb,c = 123. This 
value is similar to that obtained for fast-MC events gen- 
erated according to Model C, and it is very different from 
that of events generated according to Model B. 

We investigate the stability of our results and estimate 
systematic errors from the following studies. We divide 
the total sample into disjoint subsamples according to 
D charge, bins of "p\, xp, and Kim invariant mass, and 
repeat the analysis. We perform fit variations by chang- 
ing fixed parameters of the fit: background parameteri- 
zations, and the mass and width of the i^*(1680). We 
also repeat the analysis for samples selected with tighter 
and looser event selection criteria. These studies lead to 
the systematic errors quoted in the text and Table 0. The 
mass obtained for the n, and the mass and width obtained 



for the -?Cq(1430), are found to vary relatively little; e.g., 
m K varies in the range 770-860 MeV/c 2 . The width of 
the k, and the K7r + and NR decay fractions, are found 
to vary much more: T K ranges from 298-543 MeV/c 2 
and the decay fractions range from 28-63% and 31-5%, 
respectively. The largest NR fraction obtained (31%) re- 
mains substantially lower than that obtained without a 
k resonance. 

We have also studied the stability of our results with 
respect to the theoretical model. For example, we modi- 
fied the k Breit-Wigner to have a "running mass" term as 
proposed by Tornqvist |^] . We varied the momentum de- 
pendence of the k form factors. We introduced Gaussian 
form factors for all other resonant states. We varied the 
shape of the NR term. In all cases we obtained similar 
results for the k mass and width within errors; however, 
the details of the parameterizations affect the relative 
K7r + and NR contributions by up to a factor of two. 

Finally, we have checked whether other models without 
a scalar k provide acceptable fits. We tried a toy model 
(T) by replacing the n complex Breit-Wigner by a Breit- 
Wigner amplitude with no phase variation. This model 
converged to a similar mass and width (871 ± 10 MeV/c 2 
and 427 ± 23 MeV/c 2 , respectively) but with large decay 
fractions for this extra amplitude and for the NR am- 
plitude, reflecting strong interference. The fast-MC gave 
(Awt,c) = 60 (rms of 16) for an ensemble generated ac- 
cording to Model C, and (Awt.c) — —60 for an ensem- 
ble generated with toy model parameters. For the data, 
Awt,c = 45; i.e., the data prefers that the additional 
amplitude have a phase variation and not just a larger 
amplitude at low Kir mass. We also replaced the scalar 
k resonance by vector and tensor resonances to test the 
angular distribution. The vector resonance model (V) 
converged to mass and width values of 1103±45 MeV/c 2 
and 350 ± 93 MeV/c 2 , respectively, with a decay frac- 
tion of only 1.8% and a large NR fraction. The fast-MC 
gave (Awy,c) = 140 (rms of 23) for the ensemble gen- 
erated according to Model C, and (Awy,c) = —140 for 
the ensemble generated with vector parameters. For the 
data, Awy.c = 116; i.e., the data prefers that the addi- 
tional resonance be scalar rather than vector. We were 
not able to make the tensor model converge, the width 
being driven to large negative values. We also performed 
a variety of fits to study the NR shape jl6) in variants 
of Model A, i.e., without an additional scalar amplitude. 
We fitted the NR amplitude to polynomials, and we also 
allowed for different interfering angular distributions, but 
none of these fits were as good as that of Model C. 

In summary, we have performed a Dalitz plot analysis 
of the decay D + — > K~-k + ii + . We compared models in 
which the signal amplitude A is the coherent sum of a 
uniform NR term and Breit-Wigner Ktt resonances. Our 
best fit is obtained when we include an additional scalar 
resonance with a phase variation corresponding to that 
of a Breit-Wigner; this state subsequently accounts for 
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about half of the decay rate. The mass and width ob- 
tained are 797±19±43MeV/c 2 and 410±43±87MeV/c 2 , 
respectively. The fit mass and width of the iiTp (1430) 
depend on whether this additional Breit-Wigner is in- 
cluded or not. When not included, m ^*(!4 3 o) = 1416 ± 
27 MeV/c 2 and T K , {U30) = 250 ± 21 MeV/c 2 (statis- 
tical errors only), in agreement with PDG values [^). 
When included, m 



K* (1430) 



K* (1430) 



1459 ± 7 ± 5 MeV/c 2 and 
175 ± 12 ± 12 MeV/c 2 . Overall we conclude 
that the scalar contribution to A is not adequately de- 
scribed by the sum of a uniform non-resonant term and a 
Kq (1430) term. Including an additional scalar resonance 
in A results in a good fit to the data while the mass and 
the width of the -FTq(1430) appear higher and narrower, 
respectively, than previous reported results. 

We thank E. van Beveren and N. Tornqvist for use- 
ful discussions. We gratefully acknowledge the assis- 
tance of the staffs of Fermilab and of all the par- 
ticipating institutions. This research was supported 
by CNPq (Brazil), CONACyT (Mexico), FAPEMIG 
(Brazil), the Israeli Academy of Sciences and Humanities, 
PEDECIBA (Uruguay), the U.S. Department of Energy, 
the U.S. -Israel Binational Science Foundation, and the 
U.S. National Science Foundation. 



[1] E691 Collaboration, J.C. Anjos et al, Phys. Rev. D 48, 
56 (1993). 

[2] E687 Collaboration, P.L. Frabetti et al, Phys. Lett. B 

331, 217 (1994). 
[3] J. A. Appel, Ann. Re v. Nucl. Part. S ci. 42, 367 (1992); 

D. Summers et al, hep-ex/0009015| ; S.Amato et al, 

Nucl. Instrum. Methods A 324, 535 (1993); E.M. Aitala 

et al, Eur. Phys. J. direct C 4, 1 (1999). 
[4] D. Bartlett et al, Nucl. Instrum. Methods A 260, 55 

(1987). 



[5] D.E. Groom et al, Eur. Phys. Jour. C 15, 1 (2000). 
[6] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. 
B 308, 435(1993); CLEO Collaboration, S. Kopp et al, 
Phys. Rev. D 63, 092001 (2001). 
[7] J.M. Blatt and V.F. Weisskopf, Theoretical Nuclear 

Physics, John Wiley & Sons, New York, 1952. 
[8] E791 Collaboration, E.M. Aitala et al, Phys. Rev. Lett. 

86, 770 (2001). 
[9] N.A. Tornqvist, Z. Phys. C 68, 647 (1995). 

[10] CLEO Collaboration, D.M. Asner et al, Phys. Rev. D 
61, 012002 (2000). 

[11] There remains the question about how best to character- 
ize broad states near threshold. This subject is consid- 
ered in Ref. ^], and also recently by E. van Beveren and 
G. Rupp, Eur. Phys. J. C 22, 493 (2001). 

[12] E. van Beveren et al, Z. Phys. C 30, 615 (1986); S. Ishida 
et al, Prog. Theor. Phys. 98, 621 (1997); D. Black et 
al, Phys. Rev. D 58, 054012 (1998); J.A. Oiler et al, 
Phys. Rev. D 59 074001 (1999); M. Jamin et al, Nucl. 
Phys. B587, 331 (2000); CM. Shakin and H. Wang, 
Phys. Rev. D 63, 014019 (2001); R. Delbourgo and 
M.D. Scadron, Int. J. Mod. Phys. A 13, 657 (1998); 
M. Ishida, Prog. Theor. Phys. 101, 661 (1999); J.A. Oiler 
and E. Oset, Phys. R ev. D 60. 074023 (1999); F.E. Close 
and N.A. Tornqvist, |hep-ph/0204205| (2002). 

[13] A.V. Anisovich and A.V. Sarantsev, Phys. Lett. B 
413, 137 (1997); S.N. Cherry and M.R. Pennington, 
Nucl. Phys. A688, 823 (2001); N.A. Tornqvist and 
A.D. Polosa, in Heavy Quarks at Fixed Target, edited 
by I. Bediaga, J. Miranda, and A. Reis, Frascati Physics 
Series, Vol. XX (Laboratori Nazionali di Frascati, Roma, 
Italy, 2000), p. 385. 

[14] LASS Collaboration, D. Aston et al, Nucl. Phys. B296, 
493 (1988). 

[15] A new fit to the LASS data for the J = 1/2 S-wave 
Kit scattering amplitude in the elastic range (up to Kr\ 
threshold) yields m K *(i43o) — 1435 ± 5 MeV/c 2 and 

Fif*(i43o) — 279 ± 6 MeV/c 2 . Private communication 
from W.M. Dunwoodie for the LASS Collaboration. 
[16] I. Bediaga, C. Gobel, and R. Mendez-Galain, Phys. Rev. 
Lett. 78, 22 (1997) and Phys. Rev. D 56, 4268 (1997). 



